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Abstract 

We present results from the first lattice simulations of three dimensional non-compact 
quantum electrodynamics (QED3) with Nf four- component fermion flavors coupled to 
a weak Z2 chirally invariant four-fermi interaction. Results with Nf > 4 show that the 
scaling near the strong coupling chiral transition or sharp crossover is determined by the 
3d Gross-Neveu ultraviolet-stable renormalization group fixed point. Small deviations 
of the Nf = 4 critical exponents from the respective Gross-Neveu ones, hint at evidence 
for non-zero fermion mass generated by the gauge fields dynamics that might have 
been enhanced by the four-fermi coupling. It is also shown that the scaling region is 
suppressed at weak four-fermi couplings and large Nf values. Measurements of i) a 
monopole susceptibility which is the polarizability of the monopole configurations, and 
ii) the density of isolated monopoles, imply that for Nf > 1 and weak gauge couplings 
the monopoles do not affect the theory's confining properties, because they are shielded. 
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1 Introduction 



Spontaneous chiral symmetry breaking plays a significant role in both particle and con- 
densed matter physics. Of particular interest is the study of quantum field theories in 
which the ground state shows a sensitivity to the number of fermion flavours Nf. Three 
dimensional parity-invariant quantum electrodynamics with four-component spinors is 
such an interesting and challenging field theory with rich dynamics that resemble four- 
dimensional QCD and walking technicolor theories pQ. The gauge coupling has mass 
dimension one and thus provides the theory with a natural scale that plays a role similar 
to Aqcd in four dimensions. This implies asymptotic freedom, since for processes with 
momentum transfer k ^ the theory is effectively non-interacting. 

Non-trivial behavior may arise in the infrared limit, as suggested by an expansion in 
[2]. The theory is believed to exhibit logarithmic confinement of electric charges 
and chiral symmetry breaking when the number of fermion flavors Nf is smaller than 
a critical value Nfc- Super-renormalizability ensures that QED3 is free from ultraviolet 
divergences, thus making analytical calculations more transparent. Most analytical ap- 
proaches, mainly based on self-consistent solutions of Schwinger- Dyson equations (SDE) 
converge to values of Nj^ between three and five [3l H]. There are also SDE results 
which claim that chiral symmetry is broken for all values of Nf [5]. Also, a pertur- 
bative analysis of renormalization group flow in the large-iVj limit predicts Nfc ~ 6 
[6]. Recent progress on gauge invariant solutions of SDE has been recently reported in 
[7j. An argument based on the inequality fjR < fuv (where / is the absolute value 
of the thermodynamic free energy) that can be estimated by counting relevant degrees 
of freedom in the infrared and ultraviolet limits yields the prediction Nfc < | [8J; a 
result that was later challenged in Analytical calculations also predict that at Nfc 
the theory undergoes a conformal phase transition [TOj, which is a generalization of the 
infinite order Berezinskii-Kosterlitz-Thouless transition in two dimensions. 

Lattice simulations provided evidence that chiral symmetry is broken for Nf < 1.5 
[TT| [T^ [13], whereas Nf = 2 appeared chirally symmetric [T3] on lattices with 
physical extent up to Le^ ~ 90. The principal obstruction to a definitive answer has 
been the separation of scales in the theory, i.e. the fermion dynamical mass is at least 
an order of magnitude smaller than the natural scale [15] which is of the order of 
the momentum cut-off given by the inverse lattice spacing. In addition, large finite 
volume effects resulting from the presence of a massless photon in the spectrum, prevent 
a reliable extrapolation to the thermodynamic limit. Analytical results claimed that 
to detect chiral symmetry breaking for Nf > 1.5 lattice volumes much bigger than the 
ones currently used in numerical simulations are required [U [16]. Another possibility 
is that, based on universality arguments, the infrared limit of QED3 may be equivalent 
to the 3d Thirring model at its strong coupling ultraviolet-stable renormalization group 
fixed point, as both models are chirally invariant under the same U{2Nf) group. This 
universality argument can be valid provided in the Thirring model the interaction is 
mediated by a massless vector boson [T7j, which still needs to be checked rigorously in 
lattice simulations. So far, numerical simulations of the Thirring model predicted that 

A^/c = 6.6(1) m- 
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The existence of a critical number of flavors Nfc in QED3 can be explained by the fol- 
lowing semi-classical arguments discussed in [20] and also reviewed in |12] . In the large- 
Nf limit the photon propagator is modified by the vacuum polarisation diagram from 
to l/[k^ + ^Nfk] [H]. The super-renormalizable theory is rapidly damped in the 
ultraviolet regime and all interesting dynamics are expected in the infrared limit, where 
the dimensionless interaction strength scales as 1/A^/. Naively, one could deduce that the 
confining property of the Coulomb potential is screened by virtual fermion-antifermion 
pairs, because in the coordinate space the interaction is modified to 1/r for distances 
r S> (e^iVj)"^. However, as discussed in |20] the kinetic energy of fermion-antifermion 
pairs is positive and scales as by the uncertainty principle. In the infrared, therefore, 
both kinetic and potential terms scale as and it becomes a delicate question which 
dominates. Since when r — > the positive kinetic term must dominate the logarithmic 
Coulomb term, we deduce in this case the existence of an energy minimum at some 
non-zero r, implying the existence of stable fermion-antifermion bound states in the 
ground state. This semi-classical argument suggests a non-vanishing chiral condensate 
for Nf < Nfc- For asymptotically large r the massive fermions decouple, ceasing to 
screen the charge, and the logarithmically confining Coulomb potential is restored. For 
Nf > Nfc the theory is conformal, consisting of massless fermions interacting via a 1/r 
potential. 

Since the early 1990's, non-compact QED3 with or without extra four-fermi terms 
has attracted attention [211 122] because of potential applications to models of high Tc 
superconductivity. More recently, interest in this model has been revived by suggestions 
that QED3 with two fermion flavors may be an effective theory for the underdoped and 
non-superconducting region of the phase diagram of high- Tc superconducting cuprate 
compounds [23]. In this sense, the abstract theoretical problem of the value of Nfc 
assumes phenomenological importance. These results have also stimulated lattice simu- 
lations of QED3 with Fermi and gap anisotropies [23]. 

In this paper we present the first exploratory lattice simulation results of parity- 
invariant non-compact QED3 with massless fermions. This is achieved with the intro- 
duction of a four-fermi interaction in the QED3 action. The lagrangian for the continuum 
Euclidean field theory is given by: 

C = iJii^ - le-i^A^ + m)^i - ^{i'ii^i? + -/^.vF^,,, (1) 

with the index i implying a summation over Nf four-component fermion flavors. The 
introduction of the four-fermi interaction reduces the U{2Nf) chiral symmetry of QED3 
to a discrete Z2 symmetry: ipi — )■ •y^ip; ipi — )■ —ijj'y5. For computational purposes it is 
useful to introduce the auxiliary field a = g'^ipiipi, and the semi-bosonized lagrangian 
becomes: 

£ = ^.(^ _ ie^^A^ + a + m)^i + + -F^,F^,. (2) 

Both the non-compact ^25j and compact [26] lattice versions of this model have been used 
successfully to show that QED4 is a logarithmically trivial theory and the systematics of 
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the logarithms of triviahty follow those of the NJL model rather than those of the scalar 
theory. In the formulation of eq. ([2]) the a field acts as a chiral order parameter 
which receives a vacuum expectation value, proportional to the chiral condensate, in the 
chirally broken phase. The Dirac operator is now non-singular even with m = and its 
inversion is very fast. 

The three- and four-dimensional versions of the theory are very different from each 
other, as in 3(i the four-fermi interaction is a relevant operator. It is well-known that the 
"id Gross- Neveu model (GNM3), although non-renormalizable in weak coupling pertur- 
bation theory, becomes renormalizable in the l/Nf expansion [27J. At sufficiently strong 
couplings and large- iV/ chiral symmetry is spontaneously broken in GNM3, leading to a 
dynamically generated fermion mass S = (cr) » m. The critical coupling g"^^ at which 
the gap S/A[/y — )■ defines an ultraviolet-stable renormalization group fixed point at 
which an interacting continuum limit may be taken. As the gauge coupling is varied and 
the four-fermi couphng is fixed at some value g"^ < g"^^, then depending on the value of 
Nf the model is expected to undergo either a chiral phase transition or a sharp crossover 
from a strong coupling phase (where (ipip) 7^ 0) to a weak coupling phase where (ipip) is 
either zero or very small and possibly undetectable in current lattice simulations. Here- 
after, we will use the term "chiral transition" to denote either a chiral phase transition 
or a sharp crossover from strong to weak gauge couplings. Unlike what happens in 
Ad, here near the transition, the weak four-fermi term is expected to play a dominant 
role as compared to the ultraviolet-finite gauge interaction. Understanding the role of 
the weak four-fermi coupling in lattice simulations of non-compact QED3 is one of the 
main themes of this paper. Preliminary results were presented in [28j. Emphasis is also 
placed on the dynamics of monopoles (which in 3d are instantons) at both strong and 
weak couplings and try to understand whether they affect the confining properties of 
the model. The interplay of fermions and magnetic monopoles was recently studied in 
numerical simulations of the compact lattice version of this model [29]. The authors of 
[29] provided evidence that for Nf = 4 the monopole plasma persists even at weak gauge 
couplings. In non-compact QED3, however, the role of topological excitations may be 
different, because in this case the Dirac strings carry a non- vanishing contribution to 
the pure gauge part of the action [3D]. 

The paper is organized as follows. In Section 2, we introduce the staggered fermion 
lattice action and the monopole observables. In Section 3 we present results related to 
the strong coupling chiral transition such as estimates of the critical exponents (3^ and 
6 for Nf = 4 and show that they are close to the respective GNM3 values, implying 
that the transition is dominated by the ultraviolet-stable GNM3 fixed point. Small 
but systematic deviations of the values of Pm and 5 from the GNM3 exponents hint at 
preliminary evidence of nonzero fermion mass generated by the gauge field dynamics at 
weak gauge couplings. We also show that the scaling region is suppressed at weak four- 
fermi couplings and large Nf. Subsequently, we show that at strong gauge couplings the 
monopole density has a weak dependence on the values of the four-fermi coupling. We 
then study the behavior of the monopole susceptibility Xm (defined by Cardy in [21]) for 
Nf = 1, 12 and within the accuracy of our data we observe no diverging behavior in 
Xrn with the lattice volume, implying that monopoles and antimonopoles are shielded in 
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the continuum limit. Further evidence in favor of this scenario is provided by the density 
of isolated monopoles (positive magnetic charges that don't have any anti-monopoles in 
the nearest neighborhood) which decays faster with the inverse gauge coupling than the 
total density of positive magnetic charges. In Section 4, we summarize and discuss our 
main findings and also point to possible future extensions of this work. 



2 Lattice Model and Observables 

In this first exploratory study of non-compact QED3 with a four-fermi term, we have 
chosen the simplest Z2 chirally invariant four-fermi interaction which for practical pur- 
poses is preferable over terms with a continuous chiral symmetry, because the latter 
are not as efficiently simulated due to massless modes in the strongly cut-off theory. 
The lattice action using staggered lattice fermion fields Xi X is given by the following 
equations: 

^ = f E F,,{^)F'^{^)+Y. xix)Qix, x')x{x') + ^ E ^'(^)' 



where 



Ff,y{x) = af,{x) + a^ix + fl) - af,{x + 9) - ai,{x), (4) 
Q{x,x') = -^y^^r]^{x)[5x>^:,+f,Uxf,-5x>^x-f,Ul^^J + y^^ajx) +1715^^1. (5) 

M (x,x) 

The indices x, x' consist of three integers (xi, X2, x^) labelling the lattice sites, where 
the third direction is considered to be timelike. The symbol (x, x) denotes the set of the 
eight dual lattice sites x surrounding the direct lattice site x. Since the pure gauge action 
is unbounded from above, eq. (E]) defines the non-compact formulation of lattice QED. 
The ?7^(x) are the Kawamoto-Smit staggered fermion phases (— 1)^^^ i-^^m-i^ designed 
to ensure relativistic covariance of the Dirac equation in the continuum limit. The 
boundary conditions for the fermion fields are antiperiodic in the timelike direction and 
periodic in the spatial directions. The phase factors in the fermion bilinear are defined 
by = exp^ia^^), where a^.^ is the gauge potential. In terms of continuum quantities, 
o^xn = aey4^(x), P = f3s = where a is the physical lattice spacing. Because of the 
non-compact nature of the gauge fields, the non-fermionic part of the action is invariant 
under gauge transformations defined by the group of real numbers R. The fermionic 
part remains invariant under the smaller gauge group R/Z r^U{l). The four-fermi term 
explicitly breaks the U{Nf/2) x U{Nf/2) chiral symmetry of the lattice QED3 action to 
a discrete Z2 chiral symmetry. 

Performing simulations with massless fermions even with the reduced Z2 chiral sym- 
metry has substantial advantages, both theoretical and practical. The theory has the 
exact symmetry of the interaction terms, which forbid chiral symmetry breaking counter 
terms from appearing in its effective action. In addition, because of the large nonzero 
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vacuum expectation value of the a field at strong gauge coupling^j or its fluctuations at 
weak couplings, the model can be simulated very efficiently. Another advantage of sim- 
ulations directly in the chiral limit, is that we do not have to rely on often uncontrolled 
chiral extrapolations to measure the chiral condensate. 

The simulations were performed with the standard Hybrid Molecular Dynamics 
(HMD) R algorithm. We used conservatively small values for the HMD trajectory time- 
step dt and ensured that any 0{dt^) systematic errors are smaller than the statistical 
errors on different observables. Among the various parameters used the ones that are the 
most susceptible to algorithmic systematic errors are those near the Nf = 4 chiral tran- 
sition (/3 = 0.145) on the largest lattice volume 42^, with the weakest four-fermi coupling 
Ps = 16. By comparing the values of different observables obtained from simulations 
with these parameters and with dt = 0.00125 and dt = 0.0025 we found agreement 
within statistical errors. Therefore, we decided to use dt = 0.0025 for all 42^ and 32^ 
simulations. For the smaller lattices 16'^, 24'^ and Nf > 4 we found that dt = 0.01 and 
0.005 are small enough to suppress algorithmic systematic errors at strong coupling and 
weak gauge couplings, respectively. For Nf < 2 we used dt = 0.005 and dt = 0.0025 for 
simulations on 16'^ and 24^ lattices, respectively. 

The magnetic monopoles in the lattice model are identified following the standard 
DeGrand and Toussaint approach |32j. The plaquette angles 9^^ are written as 

Qfiu = Qfiu + 2-Ks^y{x), (6) 

where 6^1^ lie in the range (— vr, tt] and s^^i^x) is an integer that determines the flux due 
to a Dirac string passing through a plaquette. The integer number of monopole charges 
on the dual lattice sites x are then given by 

M{x) = ef,u\^f,Sux{x), (7) 

where A,^ is the lattice derivative and M G {0,±1,±2}. Since on a three-torus the 
number of monopoles is equal to the number of anti-monopoles we define the density of 
monopole charges as 

X 

We also measured the monopole susceptibility Xm introduced by Cardy [31j : 

Xm = -^J2(r'M{Q)M{r)). (9) 

r 

This observable is the polarizability of the monopole configurations and if the mag- 
netic charges are in a plasma phase, then Xm diverges implying that external magnetic 
fields are shielded. A finite Xm implies that monopoles and antimonopoles form tightly 
bound molecules. The observable Xm has been rarely used in simulations with dynamical 

"'^At large couplings, pure QED3 simulations are dramatically slowed down by the strong gauge field 
flusctuations. 
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Figure 1: (Color online) Fits of (xx) vs. /3 to eq. ([TO]) for Nf ^ 4 and /3s = 2, 4, 8, 16. 

fermions, because it is very noisy due to near cancellations of monopole-monopole and 
monopole-antimonopole contributions. With the inclusion of the four-fermi term in the 
QED3 action the algorithm became very efficient and Xm has been measured with an 
acceptable signal-to-noise ratio even at relatively strong gauge couplings. 

3 Results 

In the infinite gauge coupling limit /3 — )■ 0, it is known rigorously that chiral symmetry 
is broken [33] for values of Nf below a certain critical value. Simulations of QED3 with 
staggered fermions and /3 = have shown that the theory undergoes a second order 
phase transition at Nf ^ 8 with mean field theory exponents [20]. Therefore, as P 
increases, for Nf > Nfc there must exist a chiral symmetry restoring phase transition at 
some finite Pc- For Nf < Nfc, since the order parameter is very small in the continuum 
limit, the relic of the transition may persist as a very sharp crossover between weak 
and strong couplings with a tail of an exponentially suppressed ipip extending to weak 
gauge coupling^. For example, for Nf = 4, which is the flavor number used in the 
bulk of our simulations the SDE approaches predict that the value of the dimensionless 
chiral condensate Pi'^{^%p) is somewhere between zero and (9(10"'') [3]. In this study we 
assume that a transition takes place at some strong (pseudo-) critical gauge coupling /3c 
which depends on (ig and we use standard scaling relationships for a second order phase 
transition to extract critical exponents. We chose Nf = 4, instead of a larger Nf value, 
because as we show later in this section for large Nf values the width of the scaling 

^This may not be directly detectable in lattice simulations. 
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Table 1: Results from fits of (xx) vs. /3 to eq. (fTO|) . 





fitting range 






(/^c f^cross) 


2 


0.22-0.28 


0.313(2) 


0.96(3) 


0.09(1) 


4 


0.16-0.19 


0.212(2) 


0.96(4) 


0.052(5) 


8 


0.1325-0.16 


0.176(1) 


0.99(3) 


0.044(3) 


16 


0.1275-0.14 


0.1579(4) 


0.93 (fixed) 


0.031(3) 



region is suppressed. 

One of tfie main goals of tliis first set of lattice simulations of non-compact QED3 
witfi tfie additional four-fermi term is to understand tlie impact of tlie Gross-Neveu 
coupling on the chiral transition. As already mentioned in Section 1, in 3(i the four- 
fermi term becomes a relevant interaction (as opposed to the Ad theory where it is an 
irrelevant interaction) and therefore this term may play a significant role near the chiral 
transition even at very small values of (?^. Given that the pure 3c/ Gross-Neveu model 
with Nf = A undergoes a second order phase transition at /3sc = 0.835(1) [34], we chose 

= 2,4,8, 16 for the Nf = 4 strong gauge coupling simulations. All these values of 
Ps are in the symmetric phase of the Gross-Neveu model, implying that they cannot 
generate nonzero fermion dynamical mass on their own. The finite volume effects are 
expected to increase with (3s and in the limit (3s — )■ 00, on finite size lattices {xx) ~^ 0. 
In order to check the extent of finite volume effects, we performed simulations on 16^ 
and 24^ lattices for f3s = 2, on 32^ lattices for /^^ = 4, 8 and on 32^ and 42^ lattices for 
(3s = 16. It is clear from the data for (xx) versus (3 shown in Fig. [T]that the 24^ lattice 
size is sufficiently large to suppress finite volume effects for Ps = 2. Also a comparison 
of the values of (xx) from simulations with (3s = 16 on 32^ and 42^ lattices implies that 
32^ is large enough to suppress finite size effects for (3s = 4, 8, 16. We fitted the data for 
the different (3s values to the standard scaling relation of a second order transition order 
parameter: 

{xx) = a{(3c-(3f-. (10) 

The fitting range was varied to ensure that all the data near the transition that give stable 
values for (3m and (3c were included. For the weakest four-fermi coupling (3s = 16 the 
three parameter fits to eq. (fTOj) did not give stable values for (3m and (3c when the fitting 
range was varied. This is attributed to the fact that for such a very weak four-fermi 
coupling the scaling window is very narrow. Therefore, in order to obtain an estimate 
for Pc and the width of the scaling region we fixed Pm = 0.93 (the value of the 3d Gross- 
Neveu exponent [31]) and performed two-parameter fits with eq. f lTO|) . A reasonable 
fit (x^/dof = 1.3) was obtained from the data in the range P = [0.1275 — 0.140]. All 
the results from this analysis, namely the values for Pc, Pm, and the fitting ranges 
(where Pcross is a crossover coupling that signals the beginning of the scaling region) 
are shown in Table 1. It is deduced that the values of Pm are in good agreement (or 
slightly larger by 1 — 2 standard deviations) with the Pm = 0.93(3) of the Nf = A GNM3 
|34] . An analytical calculation based on SDK and large- A^/ approaches for QED3 with a 
U{1) chirally invariant four-fermi term [35] predicted that for Nf > Nfc = 112/37r^ the 
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magnetic critical exponent is /9m = (3 — 2a)/4a, where a = ^a/I — Nfc/Nf). This result 
reproduces the GNM3 exponent f3m, = 1 only when Nf 3> Njc- 

Our results for the critical coupling show that (3c decreases as the four-fermi coupling 
gets weaker and it is expected to saturate to an asymptotic value in the limit /3s — )■ 00. 
This asymptotic value is clearly smaller than /3c = 0.212(4) obtained from pure QED3 
simulations with Nf = 4 [12]. The discrepancy in the estimates of Pc between the two 
models is attributed to the presence of lattice discretization counter-terms in the pure 
QED3 effective action, which are forbidden by the introduction of the four-fermi term 
in the current model. It should be noted that for the same reason the inclusion of the 
four-fermi term in the non-compact QED4 lagrangian led to a larger scaling window 
in the direction of the gauge coupling than in pure non-compact QED4 [25]. Another 
interesting result shown in the last column of Table [1] is that the Gross-Neveu scaling 
region {f3c — f3cross) is suppressed as the four-fermi coupling becomes weaker. This can be 
understood by a combination of dimensional analysis and scaling arguments as follows. 
The fixed four-fermi coupling g"^ with mass dimension —1 introduces a new scale in 
the system. As the scaling region is approached by increasing /3, the magnitude of the 
dimensionless gauge coupling at the crossover into the scaling region e^j-osstt (a is the 
lattice spacing) becomes comparable to the magnitude of the dimensionless four-fermi 
coupling gla~^. Therefore, e'lj-oss'^^'i ~ qI^ where rrie is the electron mass which is an 
inverse correlation length and obeys the scaling relation: 

m^'^ = c^(/5,)(/3c-/3„oss)^ (11) 

with z/ = 1.0 for the Nf = 4 GNM3 [31]. In the large (3s limit, d{(3s) and /Across saturate 
to certain values, implying 

{(3c - /Across) ~ a/J. (12) 
V Ps 

It can be seen from the last column of Tabled] that the data for the scaling region width 
comply relatively well with eq. (fT2|) . especially for the two weakest four-fermi couplings 
= 8, 16. 

We also performed simulations on 32'^ lattices with nonzero fermion bare mass in the 
range m = 0.005, 0.020 at the respective critical gauge couplings for (3s = 2,4,8,16. 
The results for the chiral condensate as a function of m were fitted to the standard 
scaling relation: 

(XX) = cm'/'. (13) 

The data and the fitting functions are shown in Fig. [2] and the values of the critical 
exponent 6 are presented in Table [2j The extracted values of 6 are close (within 1- 
2 standard deviation) to the Nf = A GNM3 5 = 2.24(11) [31] ■ The slightly larger 
values of 6 in Table [2] as compared to the pure GNM3 value could be attributed to 
either the accuracy with which the critical couplings were measured and/or to a small 
fermion mass generated by the gauge field dynamics. Notably the 6 = 2.45(3) extracted 
from simulations with the strongest four-fermi coupling (3s = 2 hints at preliminary 
evidence in favor of mass generation by the QED3 dynamics that might have been 
enhanced by the four-fermi coupling. Future better precision simulations may further 



9 




0.005 0.01 0.015 0.02 
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Figure 2: (Color online) Fits of {xx) vs. m to eq. for Nf = 4 and /Ss = 2, 4, 8, 16. 
Table 2: Results for exponent 6 extracted from fits to eq. (fT3|) for Nf — 4. 



/3. 


6 


2 


2.45(3) 


4 


2.36(3) 


8 


2.34(3) 


16 


2.33(2) 



clarify this issue. Nevertheless, the two critical exponents Pm and 6 are sufficient to 
define the universal properties of the model at (pseudo-)criticality, since all the other 
exponents can be estimated using standard hyperscaling relations. These results clearly 
indicate that the four-fermi interaction plays a predominant role in the scaling of the 
order parameter near the chiral transition. Although the transition may be approached 
by varying the gauge coupling alone, as mentioned in Section 1 pure QED3 is super- 
renormalizable and ultraviolet-finite, and therefore the ultraviolet-stable fixed point of 
GNM3 determines the scaling properties of the transition. This situation is similar to the 
Higgs- Yukawa model where the interaction is also super-renormalizable in 3d, and 
the four-fermi interaction determines the model's universality class [36]. An analogous 
result was recently obtained in a perturbative renormalization group theory approach 
of the graphene phase diagram [37], where it was shown that in the presence of a four- 
fermi interaction the electron charge is rendered more marginally irrelevant than at the 
Gaussian fixed point. 

In Fig. [3] we present the results for the monopole density pm versus (3 for /Sg = 
2, 4, 8, 16. We see that pm decreases monotonically and smoothly with f3 with no evi- 
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Figure 3: (Color online) vs. (3 for Nf ^ 4 and /3s = 2, 4, 8, 16. 



Table 3: Results from fits of (xx) vs. /3 to eq. (fTU)) for A'^/ = 4, 8, 12 and compliance of the width of the 

scaling region to eq. ([TH) . 



Nf 


fitting range 


/3c 




R — (/3c /Across) 


75 / / /3cross(-'V/) 

^/V TV, 


4 


0.22-0.28 


0.313(2) 


0.96(3) 


0.093(8) 


0.40(9) 


8 


0.21-0.24 


0.269(3) 


0.95(6) 


0.059(5) 


0.36(7) 


12 


0.140-0.1575 


0.184(3) 


0.97(11) 


0.044(5) 


0.41(6) 



dence for any abrupt cfianges at tfie cliiral transition couplings tliat miglit signal the 
existence of a phase transition from strong to weak couplings. At strong gauge cou- 
plings pm decreases as Ps increases and above the chiral transition pm tends to become 
;9s-independent. It should be noted, however, that the impact of f3s on pm is significantly 
smaller than on the chiral condensate (see Fig. [1]); a similar result was obtained in sim- 
ulations of the compact version of the model [29j. As expected, a stronger four-fermi 
coupling produces a larger enhancement of the fermion condensate near the transition, 
whereas the monopole dynamics are mainly determined by the short distance fluctua- 
tions of the gauge fields. 

In order to study the dependence of the scaling region width on Nf we performed 
additional simulations with Nf = 8, 12 and m = 0, /3s = 2 on 24^ lattices. It is shown in 
Fig. mthat as Nf increases the values of (xx) decrease, due to the enhanced screening of 
the electromagnetic interaction by virtual fermion-antifermion pairs. The data for (xx) 
were fitted to eq. f lTO|) in the regions where stable values for (3m and (3c were obtained. 
The results for Nf = 8, 12 together with the Nf = 4 results discussed earlier, and the 
fitted curves are shown in Fig. H] and summarized in Table |3l The values of Pc decrease 
as Nf increases and the values of the magnetic critical exponent (3m are close to the 
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large- A^j GNM3 value (3^ = 1- As already mentioned in Section 1 the dimensionless 
strength of the QED3 interaction is given by ~ V^^/- Following the reasoning that led 
to the derivation of eq. f ll2p . it can be easily shown that for fixed (3s and variable Nf the 
width of the Gross-Neveu scaling region is given by: 



(/3c /Across) 



'^rossiNf) 

Nf 



(14) 



where /Across (A^/) is the crossover coupling into the GNM3 scaling region. The results 
presented in Table [3] show that the scaling region is suppressed for large Nf values. It 
is also inferred from the results in the last column of Table [3] that the data comply with 
eq. f[T^ relatively well. In a similar way to the chiral condensate, the monopole excita- 
tions density shown in Fig. [5] also gets smaller as the interaction strength is decreased 
by increasing Nf. It is also interesting to observe that in Fig. [6] the data for {xx) versus 
Pm for Nf = 4, 8, 12 and f3s = 2 collapse on a single curve, the only exception being the 
Nf = 4 data at small values of pm, a regime that coincides with the broader Nf = 4 
GNM3 scaling region. We infer from Fig. that large Nf values and strong gauge cou- 
plings produce chiral condensation and topological excitations in a very similar almost 
A^/-independent manner, without this implying anything about the impact of monopoles 
on the theory's confining properties. From the monopole density alone one cannot reach 
conclusions regarding the confining properties of the theory. 

Next, we take a closer look at the monopole dynamics by discussing results for the 
monopole susceptibility Xm defined in eq. (Q. This observable is far more informative 
than Pm, because it measures the polarizability of the monopole configurations and it 
is expected to show a diverging behavior if the monopoles and antimonopoles are in a 
plasma phase. As mentioned in Section 2, Xm is in general very noisy due to near cancel- 
lations of monopole-monopole and monopole-antimonopole contributions. However, the 
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Figure 9: (Color online) Positive monopole charges density p+ and isolated positive monopole charges 
density p\ vs. /3 for Nf = 1,4, ^5 = 2 on 32^ lattice. 

four-fermi term introduced in the QED3 action substantially increased the efficiency of 
the simulation algorithm and enabled us to measure Xm with acceptable signal-to-noise 
ratio even at relatively strong couplings. The data for Xm versus /9 obtained from simula- 
tions with N ^ = 4, = 2 on 16^, 24^ and 32^ lattices are presented in Fig. [71 Although 
Xm increases monotonically with the gauge coupling, there are no signs of a divergent 
behavior. These results imply that monopoles are shielded by forming tightly bound 
molecules. We also performed simulations with different numbers of fermion flavors in 
the range = 1, 2..., 12 on 16^ and 24^ lattices with 13^ = 2. The values of /3 = 0.4, 0.5 
used are in the weak coupling phase for all Nf values in the sense that (ipip) is consistent 
with zero and the effective order parameter decreases with the lattice volume. 

Again, the data for Xm presented in Fig. [S] show no signs of a diverging behavior as they 
are independent of the lattice volume within statistical errors. These results imply that 
for Nf > 1 the monopole dynamics do not affect the confining properties of the theory's 
continuum limit. 

The conclusion of the previous paragraph is strengthened by the data in Fig. [H where 
we show for Nf = 1 and 4 the density of positive magnetic charges (monopoles) p+ and 
the density of monopoles without any antimonopoles in their nearest neighborhood {xi ± 
1) p'+ versus beta. We fitted these data to an empirical function /(/3) = aiexp(— a2/9) 
and the extracted values of a2 and a'g for p+ and p'_^, respectively, are shown in Table HI 
The fact that p'^ decays faster than p+ confirms the scenario that at weak couplings 
the monopoles are shielded. The increase of Xm and p'_^_ at strong gauge couplings 
and/or small Nf values can be easily understood in terms of the renormalization group 
invariant Dirac quantization condition eg = engn [38] {g and gn are the monopole bare 
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Table 4: Results from fits of p+ and p'^ vs. /? to /(/3) 



aiexp(— a2/3) for — 1,4. 




1 



4 



02 
^2 



10.9(1) 
17.90(9) 



11.16(4) 
18.74(3) 



and renormalized charges). As A^/ (or /3) decreases the interaction between fermions 
and antifermions gets stronger, implying that the monopole-antimonopole attraction 
gets weaker. 

The monopole dynamics in the compact and non-compact formulations may be dif- 
ferent at weak couplings and small N f values, because in the non-compact version the 
Dirac strings carry a non-vanishing contribution to the pure gauge part of the action 
|30] . The two formulations, however, become identical in the infinite gauge coupling 
(/3 — )■ 0) and/or large- A^j (A^j — > oo) limits, as in these limits the fermionic sector with 
compact gauge links in the action dominates the dynamics. The results from simulations 
of compact QED3 coupled to a four-fermi term with N j = 4 ^9j provided evidence that 
the monopole plasma phase persists even at weak gauge couplings. Although our simu- 
lations with Nf = 4 were performed at relatively strong couplings deep in the chirally 
broken phase, it is unclear whether we reached the limit where the monopole dynamics in 
the two formulations should become similar. Furthermore, it has been suggested [39|H0] 
that the compact and non-compact formulations may be equivalent. This suggestion 
was largely based on comparisons of chirally extrapolated data for (xx) versus pm from 
the two models. The authors of [40j also suggested that even at weak couplings the 
monopoles in non-compact QED3 may exist in a plasma phase. Our results and their 
difference from those reported in [29j imply that these suggestions may not be valid. 
Marston [JT] and later on Kleinert and collaborators |l2l |33j suggested in a series of 
analytical papers on compact QED3, that above a certain critical number of fermion fla- 
vors (Nfc) a phase exists where the monopole-antimonopole potential is modified from 
1/r to ln(r), implying that monopoles are shielded and do not influence the continuum 
limit. Recently, it was reported that Nfc = 36 [13]. Also Herbut and collaborators [H] 
claimed that in compact QED3, the interaction among magnetic dipole pairs restores 
the Coulomb potential at large distances, and the monopoles are in a plasma phase at 
all /3 < 00 and Nf < 00. In the near future we plan to simulate compact QED3 with dif- 
ferent Nf values and compare the results for monopole dynamics with those mentioned 
above. 



Three dimensional QED is an interesting field theory due to its similarities to four- 
dimensional QCD-like theories and its apphcations in high temperature superconductiv- 
ity. In this paper we presented the first results from lattice simulations with massless 
fermions of non-compact QED3 coupled to a weak Z2 chirally invariant four-fermi inter- 
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action. Below we summarize and discuss our main findings. 

The values of the critical exponents (3^ and 6 extracted from simulations with Nf > 4 
are close to the values of the respective 3d Gross-Neveu exponents. This implies that 
even for very weak four-fermi couplings (almost twenty times smaller than the pure 
Gross-Neveu critical coupling) the strong coupling chiral transition/crossover for Nf > 4 
is dominated by the Gross-Neveu ultraviolet-stable renormalization group fixed point. 
The gauge interaction is an irrelevant operator near the transition, but the system can 
be driven towards (pseudo-)criticality by varying the gauge coupling. It is noted that 
certain values of the critical exponents (especially S) extracted from simulations with 
the strongest four-fermi coupling, fSg = 2, are slightly larger than those of the pure 
Gross-Neveu exponents, which hints at preliminary evidence for nonzero fermion mass 
generated by the gauge field dynamics. Future large scale simulations will further clarify 
this issue. In the weak four-fermi coupling limit the Nf = 4 critical gauge coupling /3c is 
significantly smaller than its respective pure Nf = 4 QED3 value, possibly because the 
four-fermi term forbids lattice discretization counter-terms that may exist in the pure 
QED3 effective action with m ^ 0. Furthermore, for a given Nf value, in the weak four- 
fermi coupling limit the scaling region is suppressed by a factor ~ l/y^, whereas for 
fixed (3s and variable Nf the scaling region is suppressed by a factor ~ \/ (3cross{N f) / N f . 
Also, the effects of the four-fermi coupling on the monopole density are smaller than on 
the chiral condensate, as pm depends more strongly on the short distance fluctuations 
of the gauge field, whereas the fermion condensate is significantly enhanced near the 
transition by the four-fermi coupling. For Nf > 4 we have also seen evidence that 
at strong couplings the chiral condensate is correlated to the monopole density in an 
A^y-independent manner. 

We also reached the conclusion that the monopoles are shielded at weak couplings 
for Nf > 1 because: i) the monopole polarizability measured by Xm showed no signs of 
a divergent behavior, and ii) the exponential decrease of the isolated monopoles' density 
with P is faster than the decrease of the total density of positive magnetic charges. A 
comparison of our results with those from simulations of compact QED3 coupled to a 
four-fermi term [29j which favor survival of the monopole plasma at weak couplings for 
Nf = 4 implies that the monopole dynamics in the two models are different. 

The various results from this first exploratory study are promising. This is due to the 
presence of the four-fermi coupling in the model's action, which increased substantially 
the efficiency of the simulation algorithm. In the near future we plan to extend this 
work in the following directions: i) Large scale simulations of the non-compact version 
with Nf < 4 to compare the scaling properties at the transition/crossover with those 
of pure Gross-Neveu model. Possible deviations may provide evidence for the existence 
of fermion mass generated by the gauge field dynamics, ii) Landau gauge-fixed simula- 
tions to measure the fermion dynamical mass, which according to SDE approaches it is 
expected to be an order of magnitude smaller than the natural cut-off scale in the 
continuum limit [15], and hence significantly larger than the chiral condensate, iii) Sim- 
ulations of compact QED3 with different values of Nf and comparison with the existing 
analytical results. 
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